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STEADY-bTATE UNBALANCE RESPONSE OF A THREE- DISK 


FLEXIBLE BOTOR ON FLEXIBLE, DAMPED SUPPORTS 
by Robert £. Cunnin^am 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

ABSTRACT 

Experimental data are presented for the unbalance response of a flexible, 
ball bearing supported rotor to speeds above the third lateral bending critical. 
Values of squeeze film damping coefficients obtained from measured data are 
compared to theoretical values obtained from short bearing approximation over 
a frequency range from 5000 to 31 000 cydes/min. Experimental response for 
an undamped rotor is compared to that of one having oil squeeze film dampers 
at the bearings. Unbalances applied varied from 0.62 to 15.1 gm-cm. 

INTRODUCTION 

The continuing search for improved performance of turboshaft machines is 
placing more emphasis on higher operating speeds (ref. 1). In addition to these 
higher rotating speeds there is a need to increase specific thrust ratios by re- 
ducing the overall wei^t of powerplants particularly for such applications as 
small helicopters, remotely piloted vehicles (RPV's) and small private and 
commercial aircraft (ref. 2). 

Lower specific fuel consumption is an immediate goal of engine manufac- 
turers. One area that is receiving the attention of engineering research is ex- 
cessive leakage of air between compressor stages. To reduce such losses it is 
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necessary to reduce the seal clearances, particularly those at compressor blade 
tips. These reduced clearances, however, will also mean reductions in allow- 
able amplitudes of rotor and frame vibrations. 

Lighter powerplant structures, hi^er rotor operating speeds and smaller 
blade tip clearances places a great amount of emphasis on the dynamics of 
rotors, their bearings and the manner in which the bearings are mounted in the 
frame. The combination of lighter rotors and higher operating speeds can result 
in severe vibrations when operating through one or more lateral bending critical 
speeds. 

Minimizing the amplitudes of rotor vibration at the lateral bending criticals 
means that techniques of balancing flexible lotors and shafts had to be greatly 
improved. This ultimately led to the devdopment of new and improved balancing 
procedures, one of which is the influence coefficient method as described in ref- 
erence 3. 

Sometimes, however, it is not possible to achieve the desired low levels of 
vibration by balancing alone, particularly with the built up ball bearing supports 
and rotor assemblies which comprise the modern tuxbojet engine. Rotor assem- 
blies, thou^ balanced well initially, degrade with continued service time. 

Other means had to be developed to control both steady state and transient vi- 
brations. 

Coulomb (dry friction) and viscous oil shear types of mechanical dampers 
have been used for many years to control vibration amplitudes. Normally they 
were used in controlling torsional vibrations in all types of rotating machinery 
(ref. 4) . Theinherent damping properties of hydrodynamic oil lubrication bear- 
ings were investigated by Hagg and Sankey in 1956 (ref. 5) . 

It is a well-known fact that the ball or roller bearings contain very little in- 
herent damping. Such bearings are used almost exclusively in aircraft turbine 
engines. The availability of lubricating oil and the requirement of large amounts 
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of energy dissipation in a limited space were natural inducements to the ultimate 
development and use of an oil squeeze film dampter (ref. 6) . 

To properly design such a damper it is necessary to know something about 
the dynamic characteristics of the rotating system. Normally the machinery 
design engineer will have a number of computer programs available to him. 

Such programs can provide him with the flexible critical speeds and mode shapes 
of a particular rotor model in graphical output (ref. 7) . The computer program 
referenced will also draw, to scale, a cross section of tbe rotor. It will al^ 
automatically plot critical speeds as a function of support stU^ness. Unbalance 
response computer programs are also available that are capable of predicting 
the response of almost any type of rotor on flexible, damped supports. Damp- 
ing coefficients can be an input to the program eitiier as constants or as poly- 
nomial functions of rotor speed (ref. 8). The output of such programs, however, 
are only as accurate as the values of ir.put such as support stiffness and damping 

A number of analytical papers have been written which develop the equations 
for calculating squeeze film damping and stiffness coefficients; two such papers 
which use the short bearing approximation are references 9 and 10, It is then, 
one of the objectives of this paper to experimentally determine llie oil squeeze 
film damping coefficients and to compare them with theory. 

A second objective of this investigation is to examine the experimental re- 
sponse of a three-disk rotor to varying amounts of unbalance with and without 
squeeze film damping at the bearing supports. The results then would be com- 
pared to an unbalance response computer program similar to the one described 
in reference 8. This program however, was modified first, to consider only 
circular orbits and two, an Iterative solution was added that includes support 
damping as a function of journal eccentricity ratio. 

Unbalances tested varied from a low of 0. 62 gm-cm to a high of 15. 1 
gm-cm. Tiree different types of unbalance distributions were investigated and 
these are: 



Type I - Identical set screws were placed in the same radial plane in each 
of the three disks. The set screw' in the center disk, however, 
was out of phase by 180®. 

Type II - Identical set screws were placed in the same radial plane in each 
of the two end disks only. The two screws were 180® out of 
phase. 

Type ni » Identical set screws were placed in the same plane in all three 
disks. 

Results were obtained over a speed range from 5000 to 31 000 cycles/min. 

APPARATUS 
Mechanical Features 

Experimental data was obtained with the apparatus shown in the {Aotographs 
and drawing of figures 1 throu^ 5. The apparatus consists of a long slender 
shaft supported in ball bearings which are 48. 3 cm apart. A single integral 
disk is located at midspan while removable disks of identical size, one at 
each end are located outside of the bearings (fig. 2). The rotor is 69. 6 cm 
long and has a uniform diameter of 2. 54 cm between bearings. The disks are 
9. 5 cm in diameter and contain threaded holes located every 10®, on a bolt 
circle of 8.26 cm. Set screws can be placed in the holes as correction weights 
for balancing or to deliberately unbalance the rotor in a particular bending mode. 

The rotor is supported in 204 series, radial ball bearings, one of which is 
shown in the f^otograph of figure 4. The outer races of these bearings are 
located in cylindrical steel housings, which also serve as the damper journal 
figures 3 and 4. Centering springs of 6061- T6 Aluminum shown in figure 4 
have a yield stress of 296 MN/m and when assembled locate the damper journal 
in the center of the damper bearing, this assembly can be seen in figure 3. 

The squeeze film damper used in this investigation has an (L/D) ratio of 
0.094 and is similar in design to the one described in reference 11, with the 



5 


following exceptions. No piston ring type of seals were used at the damper 
ends. The grooves, however, were retained (see fig. 3) Also given in figure 3 
are all pertinent damper dimensions. The piston ring seals were eliminated be- 
cause they introduced a tare damping and stiffness that could not be readily 
quantified. Eliminating the end seals, however necessitated a reduction in 
damper radial clearance from 127 fxm to a nominal clearance of 64 fun in order 
to keep damper oil flows within the pumping capacity of the oil supply system. 
The lubricacing and damping oil was a typical turbojet engine oil. Type II ester 
meeting the specifications of MILr-L-23699A. 

The flow of oil to each damper could be controlled independently with needle 

valves located in each of the supply lines. This flow, however, was held con- 
-4 3/ 

stant at 5. 5x10 m /min in all the tests conducted. Temperature of the oil 
could be maintained within 2° C over the elapsed time of a typical test run. 

The test rotor was driven by a 10^ 16-cm diameter air turbine having a 
maximum operating speed of 60 000 rpm. The turbine was connected to the 
test rotor by means of a 0. 95- cm diameter quill shaft and a spline coupling 
(see figs. 1 and 2) . 


Instrumentation 

Deflections of the test rotor from a neutral axis were sensed by noncon- 
tacting, eddy current probes. These probes were located in mutually perpen- 
dicular axes at each of the three disks. Distance sensor pairs were also 
cemented into the damper bearings flush with its surface to measure damper 
journal motion. 

Transducers used for measuring the forces transmitted through the damper 
oil squeeze film are quartz crystal, piezoelectric load washers. The load ca- 
pacity Is 35. 6x10^ N with a stiffness of 175x10® N/m while having a resonant 
frequency of 50 000 hertz. Linearity over the load range was 1 percent, 
these transducers can be seen in figure 5. 
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DISCUSSION OF RESULTS 

The results of this Investigation are presented in figures 6 through 14 and 
in table I, and are discussed under two separate headings. The first examines 
the response of a &ree>di8k rotor to various magnitudes of unbalance. Three 
different types of unbalance distributions were tested, these were described in 
the INTRODUCTION, they are also shown schematically in table I. The sec- 
ond part of the discussion will compare experimentally determined values of 
damping coefficients with those predicted by theory, and for unbalances ranging 
from 0. 62 to 15. 1 gm-cm and rotor frequencies to 31 000 cycles/min. These 
results are plotted in figures 13 and 14 

Unbalance Response 

Before any tests were conducted it was necessary to determine the un- 
damped critical speeds and mode shapes. These were determined from the 
computer program of reference 7 and are plotted in figure 7 for a spring sti^- 
ness of 2. 36 MN/m. The deflections are only relative but do indicate that oil 
squeeze film dampers placed at the bearing supports should be effective when 
operating through all four bending critical speeds. 

Before any deliberate unbalance distributions could be tried, the rotor had 
to be dynamically balanced. First attempts to balance the rotor without any 
squeeze film damping at the supports indicated a high degree of sensitivity of 
a ball bearing supported rotor to very small amounts of unbalance. Using a 
combination of the influence coefficient method, reference 3 and a single plane, 
vector plot, resulted in a well-balanced rotor capable of negotiating all critical 
speeds below 31 000 cycles/min with no more than a lO-jun or less, peak-fco-peak 
amplitude. 

Figures 7, 8, and 9 compare the undamped rotor response with the response 
when oil squeeze film dampers are activated. The plotted traces in figures 7,8, 
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and 9 are that of the filtered signal generated by either proximity distance sen- 
sors or piezoelectric crystal force transducers. The 10-hertz bandwidth filter 
tracks on the rotor frequency and, ttierefore, the trace is that of the rotor syn- 
chronous response only. Comparisons between undamped versus damped re- 
sponse are for very small unbalances of 0. 62 and 0. 93 gm-cm because oper- 
ating the rotor without any external damping at the bearings could have resulted 
in damaging the rotor and stationary seals if larger unbalances had been tested. 

The arr plitude versus frequency plot of figure 7 was produced by an inline, 
alternating phase distribution of weights in the three disks. The unbalance 
weights in the enr disks are 180^ out of phase with the weight in the center disk 
as depicted by the inset drawing. The total unbalance distributed over three 
disks is 0. 93 gm-cm The three set screws used as unbalance weights were 
carefully matched in weight and did not vary more than 4 mg. 

The amplitudes of motion plotted in figure 7 are of the number 3 disk. A 
maximum peak- to- peak amplitude of approximately 100 above machining 
runout, was recorded for the rotor operating with no externally applied damping. 

The undamped rotor exhibits three distinct amplitude buildups at frequen- 
cies close to the ones predicted by the critical speeds program of reference 7. 
The maximum undamped resonant amplitude occurs at about 9500 cycles/min. 

The type of unbalance distribution used in these trials would normally excite 
the first and third bending critical speeds as can be seen from the mode shapes 
of figure 6. Therefore, it is reasonable to expect larger amplitudes at resonant 
frequencies corresponding to those modes. It should be pointed out that even 
thou^ no external squeeze-film damping was present in this trial, it is obvious 
that some amount of damping is present in the oil lubricated ball bearings. The 
maximum amplitude predicted by an unbalance response program, for an assumed 
ball bearing damping of 3500 N- sec/m is 381 ^m at 7500 cycles/min. 
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Turning our attention to the lower portion of figure 7, where oil squeeze 
film dampers were activated, shows one distinct buildup »n amplitude instead of 
three. The mode shape for this damped critical occurring at 10 000 cycles/min 
was determined tx> be similar to the one at the third bending critical shown in 
figure 6. The maximum peak-to-peak amplitude, when fl:e dampers were acti- 
vated, was approximatelv 90 j^m. The peak-to-peak amplitude at the center disk 
was 216 jian, not a significant improvement for this unbalance distribution. The 
unbalance response program predicted a maximum peak-to-peak amplitude at the 
center disk of approximately 85 jim at 9600 cycles/min and about 45 (im for the 
number 3 disk. The agreement between experiment and theory for this unbalance 
distribution is not very good. Because of the large rotor excursion exhibited for 
unbalances of 0. 93 and 1. 67 gm-cm tests at greater unbalances were not attempH 
ted for fear of damaging rig components such as the labyrinth seals located at 
the bearing supports. 

The number 3 disk wants to move out at amplitudes greater than indicated 
in the undamped test of figure 7 . However, the damper journal can move out 
only as far as the clearance in the damper bearing will permit, and this is ap- 
proximately 127 pim total movement. This may account for the jagged appear- 
ance of the trace at the second and third criticals. With the activation of the 
dampers the trace smooths out; again the theoretical prediction of the maximum 
amplitude of 44 j^m is low. 

Low speed measurements made along the shaft from bearing to bearing in- 
dicated a bow or warping of the shaft was present in the ”as received^' condi- 
tion. The amount of warp or bow measured was approximately 10 to 13 ptm 
maximum at either side of the center disk. This shaft runout can be seen at 
the low speeds in all of the plots. Inadvertently, the inline unbalance distribu- 
tion was placed in a plane only 40° away from the plane of the residual warp. 

This then could account for the shaft having a strong preference to bend a greater 
amount in this plane over any other. 
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The rotor response shown in figure 8 is for an inline, alternating phase dis- 
tribution, for two disks only; the center disk did not have an unbalance wei^t 
(see inset figure). This type of unbalance distribution should have excited the 
second mode shape (see fig. 6). The undamped response, however, is similar 
to the previous response for three disk inline, alternating phase. Activating 
the oil film dampers resulted in a significant reduction in the resonant ampli- 
tude. An amplitude of approximately 12 pm over the machining runout was re- 
corded at 9700 cycles/ min while the unbalance response program predicted a 
peak-to-peak amplitude of only 1. 3 pm at 10 000 rpm. 

The unbalance distribution used to generate the experimental response 
shown in figure 9 for the number 3 disk, is an inline, inphase distribution using 
all three disks, Type III (see inset drawing). Theoil squeeze film dampers 
cause a reduction in the peak-to-peak amplitude of over 80 percent at the end 
disk and approximately 70 percent for the center disk. The damped response 
for inline, inphase unbalance distribution applied to the test rotor correlated 
well with the predicted response at very small unbalances of 0.93 gm-cm as 
shown in figure 9. 

It was decided to examine the response for much greater magnitudes of un- 
balance. In figures 10(a) throu^ (d) comparisons are made between the actual 
and theoretical damped rotor response. It is quite apparent in figure 10(a), 
which is the peak-to-peak amplitude of the center disk, that squeeze film damp- 
ers located at the bearings of a long slender rotor are going to have little effect 
on attenuating midspan excursions. The agreement between the actual response 
and the predicted is quite good for both the amplitude and the frequency at which 
resonance occurs. 

The theoretical response program will predict, with considerable accuracy, 
the amplitudes, phase an^es, and forces transmitted providing the input data, 
such as support flexibility and damping are accurate. In this case the matchup 
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at the center disk is close. Stiffnesses of the centering springs were carefully 
measured and this value (2. 36 MN/m) entered in the program as input. Damp- 
ing coefficients were calculated from theory and for this case agreed quite well 
with measured values at least at the resonant frequencv. Damping values of 
6300 N* sec/m are shown on the plots and were calculated for a concentrically 
positioned damper journal (ti = 0) . 

It should be noted here that the recorded amplitudes of motion and trans- 
mitted forces were in the vertical or y direction. ResponLic was al&o measured 
in the X-direction and appeared similar to Y measurements. The value of radial 
clearance used in calculating the damping for the theoretical unbalance response 
was 63. 5ftm, which was also the measured assembled clearance. 

Considerable disagreement is apparent between the measured amplitudes 
of motion and the theoretical for both the end disk and the damper journal, as 
shown in figures 10(b) and (c), respectively. The thereotical response is much 
greater than the measured response at frequencies above 10 000 cycles/min. It 
is not exactly clear as to why such disagreement exists at the higher unbalances 
of 10. 6 and 15. 1 gm-cm. It is apparent that the rotor end and damper journal 
opposite the turbine appear to deflect much less than predicted. The reason 
for this may be attributed to the non- cylindrical shape of the damper bearing 
(Fee fig. 5). The walls of the bearing in those locations where distance probes 
are me. .mted, are reduced in thickness. It is possible that flexing of this mem- 
ber was taking place, particularly with the higher rotating unbalance loads and 
therefore the absolute damper journal motions measured were less than if the 
damper bearing were perfectly rigid. The transmitted force also would be less 
in the actual system (see fig. 10(d)) since it is the vector sum of the elastic and 
damping film forces, both of which in turn are functions of the displacement. 

The photographs of figure 11(a) are of the damper journal orbital motion. 
The pictures on the left are for an undamped system, damper oil supply off. 
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Those pictures on the rlg^t are for identical conditions of unbalance, magnitude, 
and frequency except that now oil is being supplied to the dampers. It should be 
pointed out that these pictures are of die unfiltered signals and, therefore, 
Lissajous patterns of orbital motion are generated which can be nonsynchronous 
as well as synchronous. The scale shown to the side of each picture should be 
noted carefully since It is not the same in each picture. 

The large orbit at 9020 cycles/min is of the undamped journal motion. This 
is a synchronous critical with a peak- to- peak amplitude of 96.4 ym. The ampli- 
tude, frequency traces of figure 11(b) are for filtered signals made at the same 
time as the photographs. The first picture was snapped W'hile the rotor was 
approaching the second critical speed of 9600 rpm. The undamped trace of fig- 
ure 11(b) shows only resonances which are synchrcaious. Some of the multiple 
harmonics that occur in a lightly damped ball bearing supported rotor can have 
large amplitudes that can be quite destructive. These motions would not be seen 
at all if only the filtered trace were being observed. The Lissajous patterns at 
25 710 rpm, for the undamped rotor, figure 11(a) has one inward loop; this sig- 
nifies a multiple harmonic having a frequency twice that of the rotor speed. ITiis 
orbit is a result of the gravity vector of a horizontally supported rotor superim- 
posed on the unbalance whirl vector (ref. 12). The noteworthy point to be made 
'concerning the pictures (fig. 11(a)) is that oil-squeeze^film damping is very 
effective in attenuating the motions of both synchronous and nonsynchronous whirl. 

Oil Squeeze Film Damping 

Presented in figures 12 through 14 are values of oil squeeze film damping 
coefficients, the force transmitted through the oil film and dynamic eccentricity 
ratio (€) all plotted as functions of rotor frequency. Eccentricity ratio, as de- 
fined here, is the average of peak-^-peak journal displacements in tne X and Y 
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coordinate directions planes divided by the radial clearance in the damper. 
Table I contains a summary of average damping coefficients at different mag- 
nitudes of unbalance and types of unbalance distributions. 

The experimental values obtained in this investigation are averages cal- 
culated from amplitudes, forces and phase angles measured in both X and Y 
coordinate direction planes Damper journal distance probes had to be located 
in planes 45° from the force transducers (see fig. 5) . To obtain the correct 
phase angles the recorded angles had to be adjusted by 45°, Thisprohibly 
resulted in a slight error in calculating the damping coefficients, but due to 
the near circularity of orbits (see fig. 11(a)) the error introduced was not 
large. 

Values 01 damping coefficient were obtained by measuring the force trans- 
mitted (F^jj) through the oil film to the damper bearing, the displacement of 
the damper journal (e) , the angle <p by which the displacemait lags the ex- 
citing force (acts 180° from transmitted force) and the rotor frequency, oj. 

The derivation of the equations used in the data reduction is given in the 
appendix and the derived expressions are as follows; 


Damping coefficient 






e 


x,y 


U) 


Stiffness coefficient 




’’x, y 


e 


x,y 


Fq is the amplitude of the transmitted force. Values of stiffness coeffi- 
cient are not presented in this report. The values, as would be calculated 
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from the above expression depend on the cosine of the measured phase angles. 
The accuracy of the measured angle near 90° =*=1° (approx. ) This error is not 

serious when the sine of the angle is taken, since the sine does not vary signif- 
icantly for a 2° change when near 90°. However, the cosine does vary signifi- 
cantly for a change of 2° and, consequently, there is a large error introduced 
when calculating values of the oil film stiffness coefficient. 

Theoretical damping coefficients plotted in figures 12 throu^ 14 were cal- 
culated from the following equation: 



where p is oil viscocity; R is journal radius, L is journal land length, is 
radial clearance, C^/2; and tj is the conventional journal eccentricity ratio 
e/Cj^ defined as a static deflection divided by the radial clearance Cj^; 

T) is zero for the theoretical damping coefficients calculated in figures 12 
through 14. The eccentricity ratio e (dynamic eccentricity ratio) plotted in these 
same figures, is defined as the radius of the orbit e, divided by the radial clear- 
ance, Cj^. For T] = 0 and 180° film the theoretical damping coefficients appears 
as a straight line in figures 12 through 14. 

The above equation is derived from the Reynolds lubrication equation for oil 
and squeeze film subjected to a periodic dynamic loading. It is for the case of a 
180° film extent (cavitated bearing) using short bearing theory (ref. 10). 

Results plotted in figures 12(a) and (b) are for an inline, alternating phase un- 
balance distribution, witti a total unbalance of 0. 93 and 1. 67 gm-cm, respectively. 
Average values for the experimental damping coefficient (see table I) were 67 50 
and 6531 N- sec/m for unbalances of 0. 93 and 1. 67 gm-cm as compared to values 
of 6492 and 6213 N- sec/m for the theoretical, 180° film extent. The average 
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values calculated over the frequency range are sligiitly higjier than the theo- 
retical, however, for the most part, the agreement can be considered as good. 
The agreement between experiment and theory for the transmitted force and 
eccentricity ratio is not very good however, considerable deviation is appar- 
ent, particularly at rotor speeds above 16 000 cycles/min. The experi- 
mental data indicates that some sort of resonance is being approached between 
28 000 and 30 000 cycles/min. 

The unbalance response program predicts a resonance of approximately 
33 000 cycles/min and therefore does not reflect any significant buildup in am- 
plitude at these frequencies. This same increase in amplitude at 28 000 cycles/ 
min is also apparent in the damped trace of figure 7 . 

The results plotted in figures 13(a) and (b) are for inline, alternating phase 
unbalance distribution with the two end disks only, and designated unbalance 
Type n. The experimental values of damping coefficients have somewhat more 
scatter than the previous data, however, the values do not deviate significantly 
from the theoretical 180® cavitated film, as predicted by short bearing theory. 
Experimental values of damping coefficient appear to deviate from theory most 
at rotor frequencies near the critical which is approximately 9600 cycles/min. 
Phase angles, when approaching the critical frequency, are changing very 
rapidly and therefore, are subject to greater error than those measured at 
frequencies greater than the critical frequency. Also, rotor orbits at these fre- 
quencies, instead of being circular, are all elliptical. 

Agreement between theory and experiment for the transmitted force and the 
eccentricity ratio are quite good throu^out most of the fr^uency range. Ex- 
perimental values tend to increase at a faster rate at higher rotor speeds as In 
the previous plots. 

In figures 14(a) throu^ (d) the unbalance distribution used was Type in an 
inline, inphase distribution. Magnitudes of total unbalance tested varied from 
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0.93 to 15.1 gm-cm, however, only values of 2.25 to 15.1 gm-cm are plotted, 
l^e average experimental values listed in table I vary from 6. 8 percent below 
to 9. 2 (^rcent above the predicted values. The trend, however, indicates ex- 
perimentally determined damping coefficients, for these series of tests, are in 
good agreement with values calculated from short bearing theory and for a 180^, 
ca/itated film. Experiments also indicate that damping coefficients are not a 
strong function of the dynamic eccentricity ratio € but in general, remain 
fairly constant. Similar results are presented in reference 13, that is, squeeze 
film damping coefficients, althou^ greater in magnitude than 180° film theory 
remained fairly constant at dynamic eccentricity ratios varying from 0.03 up 
to 0. 4 ter an end fed damper with a length to diameter ratio, L/D = 0.34. 

The data plotted in figures 14(b), (c), and (d) are for total unbalances of 
4. 97, 10. 6, and 15. 5 gm-cm. It was thou^t that the resulting higher dynamic 
eccentricity ratios would result in slightly higher values of damping coeffi- 
cients. This did not occur in fact, an opposite trend was observed in that 
damping coefficient actually decreased. It was also observed that a hydro- 
dynamic instability was induced. This is referred to in die literature (refs. 

14 and 15) as a bi-stable operation. Briefly, this is where the damper journal 
can operate in a stable condition at either of two eccentricity ratios. The re- 
gion of frequency where this instability occurred is shown in figures 14(c) 
and (d) . The rapid increase in journal amplitude is also shown at approxi- 
mately 25 000 rpm in figure 10(c). A region is indicated since the bi-stable 
{^enomenon noted did not always occur at exactly the same frequency in sub- 
sequent trials. It is not the purpose of this investigation to examine nonsteady 
state behavior such as occurred, other than to note its effect in the reduction 
of the damping coefficient. An increase in the oil flow to the damper was also 
noted to occur at frequencies approaching the bi-stable region. It is postulated 
that the damper at this juncture became starved with an accompanying rupture 
of the oil film. This would then account for the decrease in damping coeffi- 
cients above the instability region. 
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SUMMARY OF RESULTS 

Experimentally determined values of squeeze film damping coefficients 
were in good agreement with theoretical values calculated from equations de- 
rived by die short bearing approximation for a 180^ film extent, and a concen- 
trically positioned damper journal (17 = 0). Values remained fairly constaiit at 
rotor frequencies to 31 000 cycles/ min. 

Central feed groove type of oil squeeze film dampers, located at the 
bearings were most effective in reducing amplitudes of motion at the rotor ends. 
They were, however, less effective in attenuating amplitudes of motion at the 
rotor center. A reduction in amplitude of over 80 percent was observed for an 
unbalance of 0. 93 gm-cm distributed inline, inj^ase over three disks. 

When oil squeeze film dampers were activated for an inline, alternating 
phase type of unbalance distribution, three distinct undamped resonant ampli- 
tudes were reduced to only one. The mode shape for tli s one corresponds to 
the third bending critical. 

The oil squeeze film dampers /ere effective in reducing the amplitudes of 
motion for non- synchronous whirl as well as synchronous. 

An unbalance response computer program was found to be reasonably accur- 
ate in predicting the response of the three-disk experimental rotor used in this 
investigation. The discrepancies that did occur apparently were due to the 
presence of a residual bend in the rotor having a shape similar to that which 
occurs at the first critical speed. The placement of weights in a line and near 
the plane of residual bend induced greater bending than would probably have oc- 
curred with a strai^t shaft. Conversely, locating the center weight 180° from 
the two end disks resulted in the rotor not bending as much as predicted by 
theory. 

Unbalances of 10 . 6 and 15. 1 gm-cm arranged inline, and inphase in all 
three disks produced an instability in the rotor referred to as bistable operation. 
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Values of dampii^ coefficient were reduced significantly at frequencies greater 
than the onset frequency. Increased flow rates of damper oil were observed that 
mi^t have resulted from the higher damper journal eccentricities and increased 
hydrodynamic pumping. 

Resonant amplitudes were observed with the undamped rotor which oo> 
curred at frequencies slightly greater than those predicted by a critical speeds, 
mode shape program. 
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APPENDIX - DERIVATION OF DATA REDUCTION EQUATIONS 

The force transmitted throu^ the oil film to the damper bearing, and mea- 
sured by the force transducers, is equal to the sum of the elastic force, the 
dissipative or damping force, and the oil film inertia force. 

Transmitted Force * Elastic Force + Damping Force + Inertia Force 

Because the volume of oil in the squeeze film is so small the inertia force 
can be ne^ected and: 




(X.y) y 


d(x.y) 

dt 


( 1 ) 


It is assumed that the journal precesses about the bearing center in an elliptical 
orbit and the motion is harmonic: 


X = e cos wt; y = e sin wt 
X y 


( 2 ) 


~ = -e^ CO sin cot; ^ = 


e„ CO cos cot 


(3) 


Let 


w,. 


dt ^ dt y 

be the amplitudes of the transmitted force in the X and y direc- 


tions, respectively; then: 


“ (^o) ^y "" (^o) V 

X y 


(4) 


Substituting equations (2) , (3) , and (4) , in (1) we obtain: 


F = K e cos cot - B e co sin cot 


(f^ cos(cot + <p^) 

X 


(5) 


and 



siiHCOt + <^y) 


F = K e sin cot + B e co cos cot * 

w ^ J <J 


( 6 ) 
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where <p and <p are the angles by which the displacement lags the exciting 
X y 

force, see sketch below: 



for y = 0; X = cot =* 0 and solving for the stiffnesj and damping coefficients 


Vx<» - wo* = H • 



Kyey(0)+ByeyO)(l) = 


Fj. sin(0 + <p ) 

y ^ 



( 8 ) 
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for X = 0; y = e^; wt = ir/2 and again solving for the stiffness and damping 
coefficients 


Vx<®> - Vx“® “ (*'ol * ’’x> 

A 

W, ”x 


»x = 




(9) 



(10) 


The transmitted forces in the X-Y coordinate planes are measured along 
with the journal displacements in orthogonal planes. The rotor frequency, cj 
is also measured and knowing these values, the dynamic film coefficients in the 
X-Y planes can be solved for. 



TABLE I. - SUMMARY OF EXPEISMENTAL AND THBORETICAL DAMPING OOEFFiaENTS 


AT VARYING TYPES AND MAGNITUDES OF UNBALANCE 


D ioper swmctry 
3.95 cm 
L-0.74 cm 
C|^ • 63. S |dn 


I - h/////////A 

rra] 


Type of oil; Type n ester meeting specU 
flcstiona of MIL-L-23699-A 
Viscosity St 38° C * 0. 033 (N-sec)/m^ 
Oil now e 5.15x10'"'* mV min 
Supidy press' 19.29x10^ N/m^ 
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- Test apparatus used in experiments on steady-state unbalance response. 
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Figure 3. - Schematic of oil squeeze film damper used in experiments 
for the steady-state response of a 3 disk rotor. 
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Figure 6, - LJndamped critical ^eeds and mode sha^s, centering spring stiff- 
ness. K-23.6aa(PN/m. 


T 


T 


T 


T 







Figure Hi. - Comparison oi rotor response for undamped versus 
damped conditions N " 3 disk, lype II unbalanced distribution; 
U “0.62gm-cm. 
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Figure 10. - Comparison of experlmentBl with theoretical dainpeii 
rotor response, total unbalance ■ 10.6 gm-xm. 







PEAK TO PEAK TRANSMITTER FORCE, N 





* 
























p TO p amp |im) 



ROTOR FREQ (CYC /min) 

Figure llfi)). - Comparison of damper journal motion at bearing no. 2 
with and without squeeze film damping. Type 1 unbalance distri* 
button. U ■ 0.93gm-cm; Cp • 63.5pm. 
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Figure 12. - Comparison between theory and experiment of 
damping coefficients, transmitted forces and eccentricity 
ratios as functions of r(^r frequency. 
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Figure 14. - Comparison between experiment and theory of 
damping coefficients, transmitted forces and eccentricity 
ratios as functions of rotm* frequen^. 
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Figure 14. ’ Continued. 
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